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Motivated by previous suggestions that three-body hard-core interactions in lower-dimensional 
ultracold Bose gases might provide a way for creation of non-Abelian anyons, the exact ground 
state of a harmonically trapped ID Bose gas with three-body hard-core interactions is constructed 
by duality mapping, starting from an iV-particle ideal gas of mixed symmetry with three-body 
nodes, which has double occupation of the lowest harmonic oscillator orbital and single occupation 
of the next N — 2 orbitals. It has some similarity to the ground state of a Tonks-Girardeau gas, but 
is more complicated. It is proved that in ID any system of N > 3 bosons with three-body hard-core 
interactions also has two-body soft-core interactions of generalized Lieb-Liniger delta function form, 
as a consequence of the topology of the configuration space of N particles in ID, i.e., wave functions 
with only three-body hard core zeroes are topologically impossible. This is in contrast with the 
case in 2D, where pure three-body hard-core interactions do exist, and are closely related to the 
fractional quantized Hall effect. The exact ground state is compared with a previously-proposed 
Pfafnan-like approximate ground state, which satisfies the three-body hard-core constraint but is 
not an exact energy eigenstate. Both the exact ground state and the Pfafnan-like approximation 
imply two-body soft-core interactions as well as three-body hard-core interactions, in accord with 
the general topological proof. 

PACS numbers: 03.75.-b,67.85.-d 



If an ultracold atomic vapor is confined in a de Broglie 
wave guide with transverse trapping so tight and tem- 
perature so low that the transverse vibrational excitation 
quantum is larger than available longitudinal zero point 
and thermal energies, the effective dynamics becomes 
one-dimensional (ID) [l], Q. 3D Feshbach resonances 
0| allow tuning to the neighborhood of ID confinement- 
induced resonances 0, 0] where the ID interaction is very 
strong, leading to strong short-range correlations, break- 
down of effective-field theories, and emergence of highly- 
correlated TV-body ground states. In the case of spinless 
or spin-polarized bosons with ID zero-range Lieb-Liniger 
(LL) [f| delta function repulsion g%§{xj — xi) with cou- 
pling constant gi — > +oo, the Tonks-Girardeau (TG) gas, 
the exact TV-body ground state was determined in 1960 
by a Fermi-Bose (FB) mapping to an ideal Fermi gas Q, 
leading to "fermionization" of many properties of this 
Bose system, as recently confirmed experimentally 0, Q ■ 

Under conditions of ultracold gas experiments, the 
ID two-body interaction Hamiltonian is usually well- 
approximated by the zero-range LL potential V% = 
52 J2i<j<k<N ^( x j ~ x k), an d the dimensionless coupling 
constant measuring its strength is 72 = mgijnt? where n 
is the ID number density. 72 <C 1 is the Gross-Pitaevskii 
regime and 72 3> 1 is the TG regime. A three-body 
zero-range interaction potential generalizing the LL- 
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interaction is V 3 = g 3 J2i<j<k<i<N S ( x o ~ x k)$( x j - x e ), 
and its dimensionless coupling constant is 73 = mg^/h 2 . 
Note that 73 is independent of the density, whereas 72 
has the density in the denominator. Hence, attainment 
of the two-body TG limit 72 3> 1 requires very low den- 
sity or very large gi, the reason why attainment of this 
limit is so difficult and the experiments 0, H| were tours 
de force. This suggests that attainment of the three-body 
hard-core limit 73 3> 1 may be easier than reaching the 
TG limit. Although three-body collisions are rare at the 
usual densities of ultracold gases, the magnitude of the 
three-body scattering length does not depend on density. 

Laughlin's wave function [9j for the ground state of a 
2D electron gas in a magnetic field in the lowest Landau 
level is exp{-^ ^=1 \ z ]\) Ui< 3<k <N( z 3 ~ z kT where v 
is an odd integer. It is an ansatz describing quasiparticlcs 
of fractional charge ejv in the fractional quantum Hall ef- 
fect. Here Zj = Xj +iyj arc complex position coordinates 
in the (x,y) plane, I 2 = h/mui c , and ui c is the cyclotron 
frequency. In the limit of infinitely strong interaction it 
is known IllT ] that the exact ground state for even 
AMs a closely related Pfaffian-like wave function = 

s t 4exp(-^ Ef =1 kf i) n£g(*J - 4) 2 nS4 - 4) 2 

where Sf4 symmetrizes over all ways of subdividing the 
N particles into two subsets of N/2 each. It describes 
particles with 2D three-body hard-core interactions, since 
^0 vanishes when Zj = z k = zg for all choices j < k < £. 
The strong similarity between the Laughlin wave func- 
tion and the ID TG ground state ipBa(xi, ■ • ■ ,xn) — 
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Xk\ in a harmonic trap 



^/h/muj the oscilla- 
tor length, suggests a possible close connection with the 
problem of three-bod y ha rd-core interactions in ID, moti- 
vating Paredes et al. [l3| to suggest a similar Pfaffian-like 
state as an ansatz for the ID ground state. It was previ- 
ously suggested [Io|, 14 1 that the fractional quantum Hall 
effect in 2D might lead to a class of quasiparticles obeying 
non-Abelian anyon statistics, with potentially important 
applications to creation of topologically protected qubits 
for quantum computation, and very recently experimen- 
tal evidence for such quasiparticles has been found fljj . 
Although no way is currently known for producing such 
quasiparticles in ID, this is further motivation for ex- 
ploring connections between three-body interactions in 
2D and those in ID, as suggested in (l 31 ] . 

Paredes et al. [13| assume that TV is even, and con- 



struct an approximate ground state by dividing the TV 
particles into two subsets of ^ each, assuming that 
each of these two is in an ^.particle TG ground state, 
taking the product of these two states, finally restor- 
ing the required bosonic symmetry over all TV! permu- 
tations by symmetrizing over all ways of choosing two 
subsets of y particles from all TV. They assumed that 
the TV bosons were trapped on a ring, thus requiring 
the ring-periodic TG ground state @. However, the 
TG ground state in a harmonic trap is much simpler 
H3 l. With that choice the ansatz of [13 is ^opar = 



[exp(- Y,f=i 



N/2 
i<j 



c jinnf<£ \xk- 



where S is the symmetrizer described above. The sim- 
ilarity with the above 2D Pfaffian-like state is evident. 



This state vanishes when xa 
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xe for all 1 < j < 



k < I < TV, because any such choice requires that at 
least two of these three coordinates, say x p and x q , lie 
in the same subset, giving a vanishing factor \x p — x q \ 
with x p = Xq. However, it is not an energy eigenstatc. 
Numerical calculations described in [l3| suggest that its 
error does not exceed 10%. In the remainder of this paper 
the exact ground state with three-body hard-core inter- 
actions will be constructed by combining a three-body 
exclusion principle with a duality mapping generalizing 
that of 0. 

Three-body exclusion principle and ideal gas with three- 
body nodes: The key to finding the exact ground state of 
a ID Bose gas with three-body hard-core interactions is 
a three-body exclusion principle and its application to 
generation of three-body nodes in an ideal gas of mixed 
symmetry. The Pauli exclusion principle requires that all 
TV-particle wave functions of a system of fermions must 
be totally antisymmetric under all permutations of its 
particle coordinates Xj , where Xj includes both the spa- 
tial position and any discrete internal quantum numbers. 
For a system of particles without internal quantum num- 
bers, one consequence is that all allowed wave functions 
vanish if two or more particles have the same spatial po- 



sition. Although a corollary of this is that allowed wave 
functions vanish if three particles are at the same point, 
this is an utterly trivial consequence. A much weaker 
requirement is that allowed wave functions must vanish 
if three or more particles are at the same point, with no 
restriction on the wave function if two are at the same 
point. This defines a "three-body exclusion principle", 
and is the starting point for the determination herein of 
the exact ground state of bosons with three-body hard- 
core interactions in ID. The Hamiltonian Hq of TV iden- 
tical particles in ID with no interparticle interactions in 
a harmonic trap consists of only the kinetic energy and 
the trap potential: 
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Its eigenstates are trapped ideal gas states, and the first 
goal here is to find ideal gas states of mixed symme- 
try satisfying the three-body exclusion principle. The 
single-particle energy eigenstates in the trap are the or- 
bitals <p n (x) = c n exp(— mui 2 x 2 /2h)H n (x/x osc ) where H n 
are the Hermite polynomials and c n is a normalization 
constant. Start from an unsymmetrized orbital product 
state ^ou with two particles in the lowest orbital and one 
in each of the remaining N — 2 lowest orbitals: 



^ou = <l>o(xi)<i>o(x2)<l>i(x 3 ) ■ ■ ■ 4>n-2{x n ) 



(2) 
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and this will remain true if its N arguments (x% , ■ ■ • , xn) 
are permuted in any way, and for all linear combinations 
of such products. Our goal is to choose linear combina- 
tions summing to a mixed-symmetry "model state" ^om 
which satisfies the three-body exclusion principle. This 
requirement is satisfied by choosing coefficients ±1 in 
such a way that $om is the sum of N terms, each of which 
is the product of the lowest orbital 4>o(xj), with j ranging 
from 1 to TV, and a Slater determinant with the remaining 
TV — 1 atoms occupying each of 4>q, 4>i, - ■ ■ , <Pn-2 once: 
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The antisymmetry of each Slater determinant ensures 
that if more than two particles are at the same point 
Xj = Xk — X£ = x, ^om w ih vanish, i.e. this (TV — 3)- 
dimensional hyper line is a three-body node of ^om- If 
and xe is in the neighborhood of x, then 
^om changes sign as xi passes through x, i.e., it is lo- 
cally but not globally antisymmetric about that point. 
Next, note that since Hq(x) — 1, each of the prefactors 
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4>o(xj) reduces to exp(— mu>x 2 /2K). Finally, elementary 
row and column operations and van der Mond's theorem 
can be applied as they were in [l2| to the derivation of 
the ground state of the harmonically trapped TG gas, to 
reduce each Slater determinant to a Bijl-Jastrow product 
of (N - 1)(N - 2)/2 factors (x k - xg). The final result is 

JV 

^om(xi, ■ ■ ■ ,x N ) = cxp(- mu>x 2 /2h) 



which vanishes if X2 — x 3 — x^ or x\ 
x\ = X2 = £4 or x\ — X2 = x 3 . 



X4 or 



JV 



x n o** ~ ( s ) 

j=i i<k<e<N,(k,e)^j 

where an irrelevant normalization factor has been 
dropped. 

Evaluation of ^qb by duality mapping: A wave func- 
tion ^ob with Bose symmetry (totally symmetric), which 
is an eigenstate of H with eigenvalue E of Eq. © 
when all N coordinates are different, and vanishes when 
any three are equal, can be found by a duality map- 
ing generalizing that used to find the TG ground state 
12j. Each locally antisymmetric factor (xk — xg) can 
be converted into an absolute value \x k — Xg\ by mul- 
tiplication by a signum factor sgn(xfe — xi) = +1( — 1), 
x k - xi > 0(< 0), yielding 

JV 

^ob(xi, ■ ■ ■ ,x N ) = cxp(- ^ mcjx 2 /2K) 



x E II \xk-x e \ .(6) 

3=1 l<k<i<N,(k,t)^j 

Since each factor sgn(xk — xi) is constant except for a 
sign change at Xk = xg, it follows that except at collision 
points Xk = xg, ^qb is still an eigenstate of Hq with 
eigenvalue Eq. It is the desired exact ground state with 
three-body hard-core interactions. It is instructive to 
write it out explicitly for N — 3 and N = 4. For N = 3 
Eq. © reduces to 

3 

^ob(xi,x 2 ,x 3 ) = exp(-^2mujx 2 j /2h) 

x(\x 2 - x 3 \ + \xi - x 3 \ + \xi - x 2 )\) , (7) 

which obviously vanishes if x\ = x 2 = ^3. For N = 4 one 
finds 

4 

^ob^i, ■ • • , xa) — exp(- mu)X 2 j/2K) 

3=1 

X(\X 2 - X 3 \\X 2 - X4\\x 3 - X 4 | 

+ \xi - I3HX1 - x 4 ||x 3 - Xi\ 
+ \xi - x 2 \\xi - X4\\x 2 - Xi\ 
+ \xt - x 2 \\xx - x 3 \\x 2 - x 3 \) , (8) 



Implicit two-body interactions: For consis- 
tency a three-body hard-core potential term 
V 3 = hm^+oogsEi^fc^JV^i ~ x k )S{Xj - x t ) 
should be added to H n to exhibit the three-body 
hard-core interaction. A subtlety that appears to have 
been missed previously is that in ID, the existence of 
three-body hard-core interactions automatically gener- 
ates two-body soft-core interactions as well. This is not 
true in 2D, but in ID it is an unavoidable consequence 
of the topology of the configuration space, and holds for 
any many-boson wave function satisfying the three-body 
hard-core constraint, whether the exact ground state 
of Eq. ^1 or an approximate ground state such as the 
Pfaffian-like ansatz of 13j , and it holds for excited states 
as well. To see this, note that the three-body hard-core 
constaint region Xj — Xk — xg is an (N — 3)-dimensional 
hyperline which is the common intersection of the three 
(jV — l)-dimensional hyperplanes Xj = x k , Xj = xi, and 
Xk = Xg, and these hyperplanes divide the configuration 
space into disjoint regions. On the other hand, in 2D 
the coordinates Zj are complex numbers in the (x, y) 
plane, the region Zj — Zk = Zg, has dimension 2jV — 6, 
and the two-body collision regions each have dimension 
2jV — 2, also hyperlines which are easily circumvented, 
and there is an exact ground state with only three- body 
hard cores, realized by a state of Pfaffian form 1(J 1 1 1 - 
Returning to the ID case, suppose, for example, that 
x \ — x 2 < x^ < x 3 < ■ ■ ■ , and we wish to realize a three- 
body hard-core interaction at x\ = x 2 = x 3 . Then the 
particle originally at x 3 must be moved to the left across 
the particle at X4, crossing the hyperplane x 3 = x±. This 
generates a finite two-body interaction at x 3 = X4, since 
local Bose symmetry about this plane requires that the 
wave function ^ satisfy ^(21, x 2 , x 3 , x 3 + e, £5, • • • ) = 
^[x\, x 2 , X3, x 3 — e,xs,---) for infinitesimal e. So long 
as d^(xi,x 2 ,x 3 , x 3 + e, X5, • • • )/de does not vanish, this 
requires that d^{x\, x 2 , x 3 , X4, 25, • • • )/dxi± change sign 
at X4 — x 3 , implying that there will be a cusp on the 
hyperplane x 3 = X4. The same proof applies, with the 
obvious changes of arguments, in the neighborhood of 
every hyperplane Xj — Xk- This is the typical behavior 
for eigenstates of a Hamiltonian containing a LL two- 
body interaction term V 2 = 92 J2i< 3 <k<N 6 ( x j ~ x k ), 
where g 2 is finite so long as "J does not vanish there, true 
only for a two-body hard-core interaction (TG limit). 
<?2 can be determined in the usual way by integrating 



Schrodinger's equation from 



X k — to Xj 



x k - 



or, perhaps more simply, by noting which Dirac delta 
function terms are generated when the kinetic energy 
operator acts on the absolute value factors in and 
noting that these must be cancelled by like delta function 
terms in a two-body interaction Hamiltonian V 2 which 
must be added to Hq to cancel them, in order that Eg^ 
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contain no delta functions. Since the value of \1/ on the 
hyperplane Xj = Xk depends on all the other coordinates, 
g 2 is not constant as in the usual LL interaction, but 
depends on all the coordinates: g 2 — g% (x\,--- ,Xpf). 
For the case N — 3, one finds after some algebraic 
reduction 

h 2 

g2(xi,X 2 ,X 3 ) = —(\xi - X 2 \ + \X! - X 3 \ + \x 2 - X 3 \) 1 . 

Am 

(9) 

As £3 approaches the hyperplane x\ = x 2 , <?2 — > +00, 
showing that the three-body hard-core interaction is 
already implied by the dependence of g 2 on x 3 , and 
the same argument applies for the other two hyper- 
planes. The complexity of the explicit expression in- 
creases rapidly with N, but g 2 is always of the form ^— 
times an expression || where the denominator D is just 
^os °f Eq. ([5]) with the exponential prefactor omitted, 
and the numerator 91 is a sum of products of factors 
\x p — x q \ of degree one less than the denominator. It ap- 
proaches +00 as any three particles approach the same 
point. 

Comparison of exact ground state and Pfaffian-like ap- 
proximation: The Pfaffian-like approximate ground state 
of Paredes et al. [13J is defined only for even N, and is 
obtained by subdividing the set of N particles into two 
subsets of y each, assuming that each of these subsets is 
in its TG ground state, taking the product of these two 
^-particle states, and restoring total Bose symmetry by 
summing over all ways of selecting from N. They 
assumed ring geometry, necessitating more complicated 
ring-periodic TG states, but here the simpler case of har- 
monic trapping is assumed. For the simplest case N = 4 
the result is 

4 

^opf(xi, ■ ■ • , £4) = exp(— muix 2 j/2K) 

x(\x 1 - x 2 \\x 3 - Xi\ + \xi - x 3 \\x 2 - x 4 | 

+ - Xi\\x 2 - x 3 \) , (10) 

which is to be compared with the exact ground state, 
Eq. ©. The Pfaffian-like state is simpler, being built 
from products of two factors \x p — x q \ as compared with 
products of three in (fTUj) , but © is exact whereas (flU)) is 
not an energy eigenstate, although it does vanish if any 
three particles are at the same point. It is easy to see 
that the Pfaffian-like state leads to soft-core LL two-body 
interactions depending on all four coordinates when the 
kinetic energy Hamiltonian acts on the absolute values, 
just as the exact ground state does, in accord with the 
previous topological proof. 

Discussion and Outlook: The exact ground state of 
a harmonically trapped system of N > 3 bosons in ID 
with three-body hard-core interactions was found by du- 
ality mapping from an ideal gas of mixed symmetry with 
three-body nodes. It was found to have two-body soft- 



core interaction cusps as well, and it was shown that this 
is an inescapable consequence of the topology of the con- 
figuration space of a system of N > 3 bosons in ID with 
three-body hard-core interactions. In view of previous 
suggestions of a possible connection between three-body 
hard-core interactions and non-Abelian anyons with po- 
tential applications to creation of topologically-protected 
qubits for quantum computation, it seems worthwhile to 
look for ways of producing three-body hard-core interac- 
tions in a ID ultracold Bose gas. This might be facilitated 
by the fact that the dimensionless coupling constant for 
three-body interactions in ID is a density-independent 
ratio 73 = mg 3 /h 2 , so higher densities could be used 
than are required to reach the two-body hard-core limit 
72 = mg 2 /nh 2 3> 1, which requires low densities n. 
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